Heavy Quarkonia Production in Relativistic d+A and A+A collisions at RHIC, as 
Q Measured by the PHENIX Experiment 

o 

CN H. Pereira Da Costa, for the PHENIX collaboration 

^ IRFU/Sphn, CEA Saclay, 

I-? F-91191, Gif-sur-Yvette, France 



(N 



This contribution presents the latest results on heavy quarkonia (J/V" and T) production in 
d + A and A + A collisions at RHIC at a center of mass energy per nucleon-nucleon collision 
I I .ysNN = 200 GeV, measured by the PHENIX experiment, as well as their implications in 

rS terms of understanding cold nuclear matter effects and the possible formation of a quark 

^ gluon plasma in central Au + Au collisions. 

5^ Heavy quarkonia have long been considered a favored probe to study the formation of a 

Quark Gluon Plasma (QGP) in heavy ion collisions for the following reasons: 1) due to their 

T— I large masses, they are hard probes produced during the first instants of the collision via hard 

^ scattering of partons (dominantly gluons, at RHIC); 2) their large binding energy makes them 

QQ harder to dissociate via interactions with the surrounding hadrons; 3) their production was 

\^ predicted to be affected by the formation of a QGP, provided that the temperature is high 

^^ enough, via a color screening mechanism similar to the Debye screening in QED^. 

^^ However their production is also affected by other mechanism that do not necessitate the 

f-^ formation of a QGP, and are therefore referred to as Cold Nuclear Matter (CNM) effects. Such 

T-H effects include: 1) the modification of the Parton Distribution Functions (PDF) in the nuclei, 

►;>. which relates to the fact that the number of partons at a given momentum fraction x inside 

a nucleon is different whether the nucleon is isolated or inside a nucleus, thus affecting the 

momentum distribution of quarkonia formed via the hard scattering of these partonsEElfl (note 

^ that a qualitatively similar effect exists at small x in the so-called gluon saturation region, 

which can be described in the Color Glass Condensate framework l^; 2) the energy loss (via 

inelastic scattering) of the projectile parton inside the target nucleus before the hard scattering 

leading to the quarkonia formation, denoted as initial state energy loss^; 3) the broadening of 

the quarkonia transverse momentum distribution due to elastic scattering of its parent partons 

inside the nucleus, denoted as Cronin effect"^ 4) the dissociation of the quarkonia bound state 

(or its precursors) by scattering off the surrounding hadrons^. These effects can be quantified 

by studying the production of heavy quarkonia in p + A (or rather d + A, at RHIC) collisions. 

Two quantities are used in this presentation to characterize changes to heavy quarkonia 

production in d + A and A + A collisions: 

• the nuclear modification factor R^a (or Raa), obtained by forming the ratio of the heavy 
quarkonia yield measured in d + A (A + A) collisions to the yield measured in p + p col- 
lisions, normalized by the number of binary {p + p) collisions, Acoi, corresponding to one 
d + A collision. This factor does not depend on the heavy quarkonia production mech- 
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anism (provided that this mechanism is a hard process), and is equal to unity if d + A 
cohisions can effectively be considered as an incoherent superposition of independent p + p 
collisions, as far as heavy quarkonia production is concerned. -R^a can be evaluated in 
bins of the heavy quarkonia rapidity (y) and transverse momentum, but also as a func- 
tion of the collision centrality, which is related to the distance between the target and the 
projectile nuclei centers (also called impact parameter, b). Central collisions correspond 
to small values of b, whereas peripheral collisions correspond to large values of b. Note 
that b cannot be measured directly, and that depending on how the centrality is actually 
measured, its correlation to b can be rather lose (especially in the d + A case). 

the central to peripheral ratio Rcp , obtained by forming the ratio of the heavy quarkonia 
yields between central collisions (for which maximum modifications are expected) and pe- 
ripheral collisions (which should be more similar to the p + p case) , properly normalized 
by the corresponding A'coi values. The advantage of -Rcp over R^a is that most experi- 
mental systematic uncertainties are canceled, since the same data set can be used for the 
numerator and the denominator. On the other hand, the quantitative interpretation of 
i?CP deviations from unity is more ambiguous due to the fact that peripheral collisions 
cannot be strictly considered identical to p + p. 
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Figure 1: Left: J/tp RdAu as a function of rapidity using 2003 d + Au minimum bias collisions. Right: J/ip Rcp 
as a function of rapidity using 2008 d + Au central (0 — 20 %) and peripheral (60 — 88 %) collisions. Solid lines 
correspond to calculations from R. Vogt based on the EKS shadowing model convoluted to different values for 

the J/?/; breakup cross-section crbioakup- 

Figure [l] (left) shows the published J/ip RdAu measured by PHENIX using RHIC 2003 d + An 
minimum bias data as a function of the J/t/j rapiditylSl. J ftp produced at negative rapidity (the 
gold going side) originate from gluons of rather large x in the gold nuclei, with x being the 
fraction of the nucleon momentum carried by the gluon. On the contrary, J/^ produced at 
positive rapidity originate from gluons of rather small x in the gold nuclei. A J/ip suppression 
is observed at positive rapidity, which is consistent with gluon shadowing (and/or saturation) in 
the small x regime inside the gold nuclei. The various curves shown on the figure correspond to 
model calculations that include a parametrization of the gluon shadowing (here EKSl^ together 
with a J/V' break-up cross-section ^breakup that varies from to 5 mb. These curves can be 
used to fit the data and extract the most probable value for crbreakup- The shape obtained using 
this most probable value describes the data reasonably well, within the large statistical and 
systematic uncertainties. 

Figure IT] (right) shows the preliminary J ftp Rqp measured by PHENIX using RHIC 2008 
d + Au data. This data set corresponds to an integrated luminosity of ~ 55 nb^^, which repre- 
sents an increase of about a factor 40 with respect to the 2003 data set, thus allowing to have 



more experimental points as a function of rapidity and use several centrality bins. The curves 
shown on the figure have been obtained using a similar approach to that of the left panel. Due 
to the reduced uncertainties and the larger number of experimental points, it appears that this 
approach (shadowing + crbreakup) is unable to properly describe the data, especially at the most 
forward rapidity, and that additional effects, such as initial state energy loss, must be accounted 
for in the calculation. Alternatively, one can parametrize the difference between the shadowing 
calculation (obtained with ^breakup = 0) and the data by a phenomenological rapidity depen- 
dent cross-section. Such a cross-section naturally increases for positive rapidity. A qualitatively 
similar trend has also been observed for lower energy ii /ijj measurements'^. 
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Figure 2: Left: J/t/) Raa as a function of centrality in 2004 An + Au collisions, for two rapidity bins. Right: J/^ 
T^AA divided by extrapolated cold nuclear matter effects based on 2008 d + Au data set and using different values 

for CTbrcakup as a function of rapidity. 

A significant ^/tp suppression has also been measured in Au + Au collisions as illustrated 
as a function of centrality (here A'part, the average number of nucleons that participate to 
Au + Au collisions at a given centrality) in figure VA (left) for two rapidity regions'i^. It appears 
that the suppression is larger at forward rapidity than at mid-rapidity. It is important to 
quantify which fraction of this suppression is due to CNM effects alone, in order to single out 
the effect of the possible formation of a QGP. This is achieved by using a parametrization of 
the observed d + Au suppression (using a shadowing model and a rapidity dependent iTbrcakup) 
and extrapolating it to the Au -|- Au case. The measured Raa can then be divided by this 
CNM-only expected suppression, as shown in figure^ (right )'^. A significant i/ip suppression 
remains for central enough collisions (large values of A'part) even after CNM effects have been 
removed. This indicates that additional effects (for instance due to the formation of a QGP) 
must be accounted for to be able to describe the J/f/' suppression observed in central Au -|- Au 
collisions. Besides, this remaining suppression is now similar between the two rapidity domains, 
meaning that the differences observed in figure ^ (left) can be entirely accounted for by CNM 
effects. Also shown on the figure are the results obtained in Cu -|- Cu collisions^ and for which 
no significant suppression beyond CNM effects is observed. 

Finally, T resonances have been measured by PHENIX in both p + p (using 2006 data set) 
and d -|- Au collisions (using the 2008 data set), allowing one to form the first T R^a at this 
energy. Figure p^ (left) shows the di-muon invariant mass distribution at the T mass, measured 
in one of PHENIX muon arms in d -|- Au collisions. The various coloured bands correspond to 
the different contributions to this spectra, and are fitted to the data in order to evaluate the 
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Figure 3: Left: di-muon invariant mass distribution at tlie T mass measured using 2008 d + An collisions. Coloured 
area correspond to the different contribution to the di-muon spectra. Right: T ^R^a as a function of rapidity. 



number of produced T resonances. Figure ^ (right) shows the T i?^A obtained by combining 
these data to yields measured in p + p colhsions. A suppression is observed at forward rapidity, 
similarly to the J/'0 case, which has yet to be compared to theoretical calculations. 

T resonances have also been measured in Au + Au collisions at mid-rapidity by the PHENIX 
experiment. However, due to the limited statistics, only a 90 % Confidence Level (CL) could be 
evaluated so far for .Raa and the contributions to the di-lepton invariant mass spectrum from 
other (background) sources in the T mass region have not been accounted for. The resulting 
value is -Raa < 0.64 at 90 % CL for all di-leptons produced at mid-rapidity (|y| < 0.35) with 
mass Me [8.5,11.5] GeV/c^. 

To summarize, a significant J/ip suppression has been observed by PHENIX in both 2003 
and 2008 d + Au data sets at positive rapidity (in the gluon shadowing/saturation region for the 
gold nuclei), with an increased statistics of about a factor 40 for the 2008 data. This suppression 
can be attributed to CNM effects, but cannot be properly described (for the 2008 data) with a 
simple model that considers only PDF modifications and a (rapidity independent) J/^ break-up 
cross section. Phenomenologically parametrized CNM effects can be extrapolated to Au + Au 
collisions and compared to the J ftp suppression measured in these conditions. It appears that an 
additional suppression (not attributed to CNM effects) remains for central Au + Au collisions, 
which is largely rapidity independent. Concerning the T resonance, a suppression similar to that 
of the ,l/ip has also been observed for the first time in d + Au collisions and a 90 % confidence 
level for high mass di-lepton pairs (dominantly T mesons) Raa has measured. 
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